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' Abstract: We propose a new model of holographic dark energy with an action. It is 

the first time that one can derive a HDE model from the action principle. The puzzles of 
causality and circular logic about HDE have been completely solved in this model. The 
evolution of the universe only depends on the present state of the universe, clearly showing 
that it obeys the law of causality. Furthermore, the use of future event horizon as a present 
cut-off is not an input but automatically follows from equations of motion. Interestingly, 
this new model is very similar to the initial one of Li except a new term which may be 
explained as dark radiation. 
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1 Introduction 

The dark energy problem is a longstanding problem ever since the discovery of the accel- 
erating expansion of the universe [1]. For a recent review of dark energy, please refer to 
[2]. Based on an idea of Cohen et al [3, 4], Li proposed the first model of the holographic 
dark energy (HDE) which can drive the accelerating expansion [5]. Though it is in good 
agreement with observational data [6] , it causes some criticisms due to its use of the future 
event horizon as a present cut-off. These criticisms can be summarized as the causality 
problem and the circular logic problem. The causality problem: it seems that the evolu- 
tion of the universe depends on the future information about the universe, the future event 
horizon. Beside, the equations of motion are non-local since the future event horizon is 
defined globally. The circular logic problem: the future event horizon exists only in an 
accelerating universe. How can one use an assumption based on the accelerating expansion 
to explain the accelerating expansion? For recent interesting discussions of these problems, 
please see [7]. For proposals using other infrared cut-off, see [8]- [10]. 

In this paper, we solve the causality problem and circular logic problem of HDE com- 
pletely. First, we derive a new HDE model from the action principle. We find that we can 
rewrite all equations of motion (both the new and the initial HDE models) in local forms. 
Then the evolution of the universe only depends on the present initial conditions, which 
clearly obeys the law of causality. What is more, we do not need to assume the future 
event horizon as the infrared cut-off but only a local equation L = — ^. Magically, the 
equations of motion force the cut-off to be exactly the future event horizon aL. So it is 
not that the present cut-off depends on the future event horizon but conversely the future 
event horizon is completely determined by present cut-off through equations of motion. 
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This paper is arranged as follows. In Sec. 2, we derive a new HDE model from the action 
principle. We find exact solutions with matter and radiation and investigate cosmology 
in this model. In Sec. 3, we study models with other infrared cut-off in the action. We 
conclude in Sec. 4. 



2 New HDE model from action principle 

In this section, we derive a new model of holographic dark energy from the action principle. 
We find exact solutions of this model and prove that the cut-off is exactly the future event 
horizon. We also investigate cosmology of this model and find that radiation is dominant 
in early time, while dark energy is dominant in late time. 

2.1 The general theory 

Consider the Robertson- Walker metric 

ds'^ = -N\t)dt^ + a\t)[— ^ + r^dO^] (2.I) 

1 — Kr^ 

and the following action 

« = ^5^) - mm + ^)i + sm, (2.2) 

where R is the Ricci scalar, \/—g = Na^ (we have integrated the r, 9, (p parts) and Sm 
denotes the action of all matter fields (we use M to denote all matter fields, m to denote 
matter without pressure and r to denote radiation). In following derivations, we first take 
the variations oi N,a,\,L and then redefine Ndt as dt. We obtain 



d 2 1"^ c A Svr 

2da + d'^ + K c A 

and 



3a2L2 6a4 



Sttpm, (2.3) 



L = — , L= —— + L{oo) 



a J I a{t' 

We want to mention that we shall prove L{oo) = using the asymptotic solutions to be 
derived below. So, quite interestingly, aL is exactly the future event horizon. Besides, from 
eq.(2.3), it is easy to see that the A(0) term behaves the same way as radiation, thus can 
be a candidate for dark radiation [11]. For the purpose of solving equations of motion, we 
can always redefine X{t) as A(t) — A(0), pr as pr + iQ^^i 1 to let A(0) = 0. From eq.(2.3), we 
can derive 

a A 47r 

- = - ^(^A/ + 3m , 2.5) 
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Note that we have A < for large enough time from eq.(2.4) (c > 0). We can ignore the 
matter effects for large enough time when dark energy is dominant. So eq.(2.5) implies the 
accelerating expansion of the universe. 

Let us now try to find out solutions of the above equations. Fortunately, we can solve 
eqs. (2.3-2.5) exactly with general c, k, pm oc ^ and As a simple example, it is easy 

to observe that there is a de Sitter solution with (a = e-^*, L = e~^*/H, A = -6H'^e~^^*) 
and (c = 6,K = 0,pm = Pm = 0). We postpone discussing the general case to the next 
section. 

Let us go on to discuss the asymptotic solutions for large enough time. Using eq.(2.3), 
we can derive 

d(da) c 47r 9, , , , 

^^ = ^-'^ + Ya'(m-3pAf). (2.6) 

Note that is the only increasing function of time ( L is a decreasing function) on the 
right hand of eq.(2.6), which becomes dominant for large enough time. So, we have 



d{aa) c 
dt " 3l2 

asymptotically. Applying L = — we can rewrite the above equation as 

d'^a ca 



(2.7) 



(2. 



with the general solution 
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a = ciL 2 + C2L 2 ~L 2 . (2. 



Using L = — i and the asymptotic solution a ^ L 2 ^ we can derive 



a~t='-\/^, c>6 (2.10) 

a~e^^*,C3>0, c = 6 (2.11) 



(c4-t)'-^^, c<6. (2.12) 

Eq.(2.12) implies that time would end at t = C4. Using eq.(2.9) with c < 6, we can derive 
aL — > when L — ?• or equivalently f — >■ C4. Thus the future event horizon shrinks to a 
point (aL — t- 0) at t = C4. Then time would end. In fact, below we shall see that we have 
w < —1 when c < 6. So there would be a big rip for the universe with c < 6. 
From the state parameter of dark energy 

Pde AL2-2ca2 -3 + 2C + V9 + 12C 

w = = — — K K = , (2.13) 

Pde 3AL2 + 6ca2 3(-3 - 2c + V9 + 12c) 
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we get 



w = —1, c = 6, 
w < c> 6, 

u; < -1, c<6. (2.14) 



Note that 

Pde oc a-3(i+-) > pM, (2.15) 

consistent with our above assumption that dark energy is dominant for large enough time. 

The asymptotic solutions eqs.(2.10,2.11) imply a(oo) — )■ oo and eq.(2.12) implies a{t = 
C4) — )• 00. So we have L(oo) — )■ or L{t = C4) — )• from eq.(2.9). Thus, it is very 
interesting that aL is exactly the future event horizon. It should be mentioned that the 
above method can also apply to the initial HDE model. Using asymptotic solutions, we 
need not to assume that aL is exactly the future event horizon but only L = — i. Then 
all the equations of the initial HDE model become local and depend on only the present 
initial conditions. 

Let us comment on our results. First, it is the first model of holographic dark energy 
with action principle. Second, it implies the accelerating expansion of the universe. Third, 
eqs. (2.3,2.4) show that the evolution of the universe only depends on the present initial 
conditions a, d, L, A. So, it is clear that our HDE model obeys the law of causality. Fur- 
thermore, the use of the future event horizon as a cut-off is not an input in our model. 
Instead, the cut-off aL turns out to be the future event horizon automatically from equa- 
tions of motion. So the long-standing problem of HDE "Why does the present evolution 
of the universe depend on the future of universe" has been solved. In fact, the evolution of 
the universe only depends on present conditions but equations of motion force the future 
event horizon to be the present cut-off aL magically. In a word, in this model we answer 
the above question by "It is natural that the future of the universe (future event horizon) 
depends on the present conditions (aL) of the universe." Fourth, this model includes a de 
Sitter solution. Fifth, the energy density of dark energy 

will remain positive if it is positive at the beginning and a > 0. Finally, from eq.(2.2) it is 
interesting to note that (—A) is the conjugate momentum of the cut-off L. Or equivalently, 
A is the energy associated with the conformal time t] = —L. 



2.2 The exact solutions 



In this section, we shall solve equations eqs. (2.3-2.5) exactly with matter (pm = j^jPm = 
0) and radiation (p^ = -^,Pr = ^Pr)- Thus, we have pu = + -^^Pm = g:^- For 
simplicity, we focus on the case with k = 0. It should be stressed that one can solve these 
equations exactly with general k. 
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Redefine a new time dL = —drj = we can rewrite eq.(2.6) as 

d'^a ca , , , 

Note that the above equation is independent of the parameter d. It seems that radiation 
does not affect the evolution of the scale factor a. This is however not the case. As we have 
mentioned before, we can always redefine A as (A — A(0)) , d as (d + ^^), to let A(0) = 
(note A(0) means A(t = 0) instead of A(L = 0)). Then, using A(0) = we can relate the 
radiation parameter d to the scale factor a. 
The general solution of eq.(2.17) is 

a = ciy + 02^2 + iftL^lnL, c = 6, (2.18) 
L 3 

a=ciL 2 +C2L 2 L^, c/6. (2.19 

c — 6 

Let us first discuss the case with c = 6. Using dt = —adL, we get 

i = -ci In L - -C2L^ - -bL^ In L + -^L^. (2.20) 
3 9 2 T 

From eq.(2.5), we can derive 

A = -6[a^d + ab + d] = -6[aa" - a'^ + ab + d] 

= -Qd H 45L C2 + 12L C2 jj 

-U^L^lnL + 8bL^C2lnL-^^^^^ + -b'^L^lnlf, (2.21) 
3 L 3 



with a' = ^f. From the above equation, we can derive A' while we can also derive A' = 



2^C 

from eq.(2.4). One can check that these two results are equal to each other, which implies 
that we have got the self-consistent solutions to eqs.(2.3-2.5). With eqs.(2.18,2.21), we can 
obtain the asymptotic state parameter as 

w -1, L — ;> 0, 

L^Lo, (2.22) 

where Lq is defined by a(-Lo) = 0, denoting the beginning time of the universe t = 0. 

Now let us turn to discussing the case with 6 briefly. With eq.(2.19), we can derive 

A = -6[a^S + ab + d] = -6[aa" - a'^ + ab + d] = -6d + ... (2.23) 

One can check that eqs. (2. 19,2.23) satisfy all of eqs.(2.3-2.5). From eq.(2.4), it is easy to 
observe that A(0) = limAjc^o- Applying the condition 

limA|c^o = -6(26ci + d-cf) = 0, (2.24) 

we can derive 

4 = d + 2cib, (2.25) 
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which shows that radiation does affect the evolution of the scale factor a eq.(2.19). 
One can derive the asymptotic state parameter as 

-3 + 2c + V9 + 12c 1 ^ ^ 

-1 < w — , < --, L 0, 

3(-3 - 2c + V9 + 12c) 3 

w^-^, L^Lo, (2.26) 

for c > 6 and similarly 

-3 + 2c + V9 + 12c ^ ^ ^ 

w-)-— < -1, 

3(-3-2c + V9 + 12c) 

w^-^, L^Lo (2.27) 

for < c < 6. 
2.3 Cosmology 

In this section, we briefly discuss cosmology in our model. It is necessary to check that 
radiation is dominant in early time and a ~ y/i thus consistent with the standard cosmology. 
In early time, we have 

a(Lo) = 0, a{L)f^a'{Lo){L-Lo), a'(Lo) < 0. (2.28) 

Applying L = - i, we get 

t^-^{L-Lo)'. (2.29) 
Thus, we obtain the expected result 



a{t) w v^-2a'(Lo)t , (2.30) 

implying that radiation is dominant in early time. Using the above solutions and A(0) = 0, 
we can derive 

Wde^-^, f^de~0, ~ 1, (2.31) 

which is consistent with eqs. (2.22,2.26,2.27) in the early time limit L — t- Lq. Using the 
exact solutions eqs. (2. 18, 2. 19), we can obtain the same results as above. So our model has 
passed the check that radiation is dominant in early time and a ~ y/i. Thus, it would not 
ruin standard results such as nuclear genesis. 

Applying the asymptotic or exact solutions, we can easily find that dark energy is 
dominant in the late time and state parameter behaves as 



w = —1, c = 6, 
1 

w<--, c> 6, 

w<-l, c<6. (2.32) 
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When c = 6, the universe will approach de Sitter space asymptotically. While for c < 6, 
dark energy will behave as phantom in late time and end up with a big rip. So it is 
necessary to estimate the value of c. 

Since the universe has turned to the phase of accelerating expansion just recently, now 
we have f ~ which leads to ^ ~ —^Pm from eq.(2.5). Using aL ~ together with 
Pde ~ |pmi we have ^J^i ~ ^y- ~ It follows that c 2.6 which implies that, 

similar to the initial HDE model, dark energy will behave as phantom in late time. It is 
interesting to check our model with the observational data. We leave it to the future works. 

Finally, we want to talk about the radiation-like term A(0) in our model. We can 
separate the energy density {-^^ + g^^*,^, ) and pressure (^g^ - ^A-Jii^L^ ) i^^^^ eq.(2.3) 
into the dark energy part and dark radiation part: 

A - A(0) c 

Pde = —TT. 2— + 



levra^ 87ra2L2 ' 
_ A - A(0) _ c 
~ 48^o4 " 2Wl2 ' 
A(0) 1 

Pdr — -YQ^^A ' Pdr — -^Pdr- {1.66) 

By dark radiation, we mean that it has no interactions with other fields. It originates from 
the initial condition of A and never decays. Interestingly, observational evidences support 
the existence of dark radiation [11]. We shall study this topic in details in the future work. 

3 Models with other cutoff 

In this section, we use the particle horizon and the Hubble horizon as the cut-off in the 
action. We find that, similar to the initial HDE model, neither of them can drive the 
accelerating expansion of the universe. 

3.1 The particle horizon as the cutoff 

In this subsection, we shall use the particle horizon as a cut-off. For simplicity, we focus 
on the case k = pM = Pm = 0. We find that this model has no accelerating solutions. Let 
us start with the action 

Following the same procedure of the above section, we obtain 

6 



,a 2 c A 



a 



3a-^L^ 6a4 ' 

(3.2) 



2aa + a? c 



2 3a2L2 ' 



and 



f 1 . dt' 
L = L= -— + L(0) 



a' Jq a{t' 

4ac , f'Aa{t')c 
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From eq.(3.2), we can derive 



a A , 

If c > 0, eq.(3.3) implies A < so that f < for large enough time. On the other 
hand, eq.(3.2) implies A < if c < 0. So this model has no accelerating solutions in either 
case. 

3.2 The Hubble horizon as the cutoff 

Now let us turn to the case of the Hubble horizon as the cut-off. We focus on the case 
K = 0. We begin with the following action 

Following a similar procedure as the above sections, we obtain 

a 2 c A Svr c a o Svr 

^a^ = 3l^ - 6^ + YP'' = S^-J + T^^^' ^'-'^ 

2da + c A A(d^ — aa) ^ 

a? Qaa? 'ia?cfi 

c Tda + „ „x 
= 2 Sttpm, (3.7) 



and 



3 



.^^lA.i^. ,3.S) 

From eqs.(3.6,3.7), we can derive 

d cd A-K , ^ . , ^, 

- = ^(PM + 3m)- (3.9) 

a 3 a 3 

From eq.(3.6), it is clear that (1 + |) > 0. Then eq.(3.9) leads to ^ < 0. So there are no 
solutions of accelerating universe in this model. 



4 Conclusions 

In this paper, we propose a new HDE model with the action principle. This is the first time 
that one can derive a HDE model from an action. Furthermore, the causality problem has 
been solved completely in this and the initial HDE models. By introducing two fields A and 
L, all equations of motion become local and the evolution of the universe is determined 
completely by present initial conditions. So clearly the HDE models obey the law of 
causality. Quite interestingly, the circular logic problem in HDE can also be solved. One 
always criticizes the original HDE model since the use of the future event horizon as a 
cut-off has assumed the accelerating expansion of the universe. Then how can one use an 
assumption based on the accelerating expansion to explain the accelerating expansion? In 
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fact, we do not need to assume the future event horizon as the cut-off but only a local 
equation L = — i. Then the equations of motion magically force the cut-off to be exactly 
the future event horizon. This is also a support of the causality in HDE models. As we 
have shown, it is not that the present evolution of the universe depends on the future event 
horizon but that the future event horizon is determined by the present cut-off through 
equations of motion. 

It is interesting to note that this new model is very similar to the initial one. For 
example, in both models, only the future event horizon rather than the particle or Hubble 
horizons can drive the accelerating expansion of the universe. Beside, both models have a 
pure de Sitter solution and behaves as phantom below the critical parameter( c = 6 in the 
new model and = 1 in the initial one). The main difference between the two models is 
that the new one predicts the existence of dark radiation. We shall study this interesting 
problem in the following works. 
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